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ABSTRACT 



Or 



' Aims. We present N-band spectro-interferometric observations of the red supergiant WOH G64 in the Large Magellanic Cloud (LMC) 

using MIDI at the Very Large Telescope Interferometer (VLTI). While the very high luminosity (~ 5 x 10 5 Lq) previously estimated 
for WOH G64 suggests that it is a very massive star with an initial mass of ~40 Mq, its low effective temperature (~3200 K) is in 
serious disagreement with the current stellar evolution theory. 

Methods. WOH G64 was observed with VLTI/MIDI using the UT2-UT3 and UT3-UT4 baseline configurations. 
Results. The dust envelope around WOH G64 has been spatially resolved with a baseline of ~60 m — the first MIDI observations 
to resolve an individual stellar source in an extragalactic system. The observed Af-band visibilities show a slight decrease from 8 to 
5_( . ~10 yum and a gradual increase longward of ~10 yum, reflecting the 10 /im silicate feature in self-absorption. This translates into a 

steep increase of the uniform-disk diameter from 8 to 10 fim (from 18 to 26 mas) and a roughly constant diameter above 10 yum. The 
visibilities measured at four position angles differing by ~60° but at approximately the same baseline length (~60 m) do not show a 
noticeable difference, suggesting that the object appears nearly centrosymmetric. The observed N-band visibilities and spectral energy 
distribution can be reproduced by an optically and geometrically thick silicate torus model viewed close to pole-on. The luminosity 
of the central star is derived to be ~ 2.8 x 10 5 Lq, which is by a factor of 2 lower than the previous estimates based on spherical 
models. We also identify the H 2 absorption features at 2.7 and 6 fim in the spectra obtained with the Infrared Space Observatory and 
the Spitzer Space Telescope. The 2.7 /jm feature originates in the photosphere and/or the extended molecular layers, while the 6 j-im 
• feature is likely to be of circumstellar origin. 

Conclusions. The lower luminosity newly derived from our MIDI observations and two-dimensional modeling brings the location of 
WOH G64 on the H-R diagram in much better agreement with theoretical evolutionary tracks for a 25 Mq star. However, the effective 
temperature is still somewhat too cool compared to the theory. The low effective temperature of WOH G64 places it very close to or 
even beyond the Hayashi limit, which implies that this object may be experiencing unstable, violent mass loss. 
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1 . Introduction vantage of studying RSGs in the LMC and SMC is that we can 

probe possible metallicity effects on the mass loss, because the 

The evolution of massive (M > 8 M ) red supergiants (RSGs) is me tallicities of the LMC and SMC are systematically lower than 

not well understood, because of several physical processes dif- ±e solar value ( „ .5Z© and 0.2 Z Q for the LMC and SMC, 

ficult to theoretically formulate, s uch as mass loss, convective respectively). Although radiation pressure on dust grains is of- 

mixing, and rotation (e.g., Massey |2003]l. Comparison of theo- ten considered to be the driving mechanism of mass loss in cool 

retical evolutionary tracks for RSGs with observational data is evolved stars, it is little understood where and how dust forms in 

crucial for testing the present theory of the evolution of mas- RSGs and how mass outflows are initiated. Since dust formation 

sive stars. While it is usually difficult to reliably derive the lu- is expected to depend on metallicity, studies on the circumstellar 

minosities of Galactic RSGs, their counterparts in the Large and environment of RSGs in the LMC and SMC are useful for clar- 

Small Magellanic Clouds (LMC and SMC, respectively) provide ifying the me tallicity effect on the dust formation and obtaining 

us with an excellent opportunity to compare theory and observa- insig hts into the driving mechanism of mass outflows, 
tion on the H-R diagram, thanks to their known distances (~50 

and -60 kpc for the LMC and SMC, respectively). Another ad- WOH G64 (IRAS04553-6825) is a highly lumino us coo l 

star in the LMC first identified by Westerlund et al. d!981|l. 



Send offprint requests to: K. Ohnaka Subsequent infrared observations by Elias et al. dT986j and 
* Based on observations made with the Very Large Telescope Roche et aL ^3 revealed that it exhibits a huge infrared ex- 
Interferometer of the European Southern Observatory. Program ID: cess and the 10 ^ m silicate feature in self-absorption, which 
076.D-0253, 080.D-0222 means that the object is surrounded by an optically thick dust 
** This work is based [in part] on observations made with the Spitzer envelope, experiencing heavy mass loss. WOH G64 i s also 
Space Telescope, which is operated by the Jet Propulsion Laboratory, know n to show OH, SiO, and H 2 Q mas ers (Wood et al. 1986 
California Institute of Technology under a contract with NASA. 119921 van Loon et al. 119961 119981 120011 Marshall et al. 120041 1. 
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The spectroscopic studies on the TiO absorption bands suggest 
spectr al types of M5-M7 (Elias et al. 119861 Van Loon et al. 
120051 1. which corresponds to effective temperatures of 3200- 
3400 K (see Sect. E). Elias et al. d 19861 ) and van Loon et al. 
(2005 ) derived bolometric magnitudes of -9.7 (6 x 10 5 Lq) and 
-9.4 (5 x 10 s Lq), respectively, assuming spherical dust shells. 
Figure Q] shows an H-R diagram with the evolutionary models 
for Z = 0.008 of Schaerer et al. ( 1993 ). For comparison, we also 
plot the evolutionary track for a 25 Mq star with Z = 0.01 newly 
computed with another evolutionary code (author: T. Driebe), 
which is based on the program of Blocker ( 1995 ). The convec- 
tive overshoot is included as described in Herwig et al. (1997 ) 
and Herwig (2000 ) with an overshoot parameter of / = 0.004, 
and the mass loss is incorporated by using the formula of de 
Jager et al. (1988 ) from the zero-age main sequence. Figure Q] 
suggests that a star with an initial mass of ~40 Mq can reach 
the high luminosities observationally derived for WOH G64. 
However, the evolutionary models predict that a 40 Mq star 
evolves only down to ~6300 K and does not reach the low ef- 
fective temperatures of 3200-3400 K observed for WOH G64. 
Therefore, the location of WOH G64 on the H-R diagram is in 
serious disagreement with theory. 

It is worth noting here that Roche et al. (1 993 ) and van Loon 
et al. ( 119991 ) suggest a disk-like or clumpy structure for the dust 
envelope of WOH G64. Van Loon et al. d 19991 ) attempted to re- 
produce the observed spectral energy distribution (SED) with a 
spherical dust shell model, but could not find a reasonable fit. 
The reason is that the huge mid- and far-infrared excess indi- 
cates a large optical depth for the dust shell, which would lead 
to very low flux in the near-infrared and in the visible. However, 
the observed SED of WOH G64 shows high fluxes at these wave- 
lengths. The possible presence of a disk-like or clumpy structure 
makes the above luminosity estimates assuming spherical dust 
shells uncertain. If there is indeed such a deviation from spheri- 
cal symmetry, the luminosity of the central star can only be esti- 
mated from multi-dimensional radiative transfer modeling using 
as many observational constraints as possible. Therefore, obser- 
vational studies on the disk-like/clumpy structure in WOH G64 
is important both for deriving the luminosity of the central star 
and for a better understanding of the mass loss in RSGs in the 
LMC. 

High spatial resolution observations in the mid-infrared are 
a powerful tool to obtain a clearer picture of the circumstellar 
environment of RSGs. The MID-infrared Interferometric instru- 
ment (MIDI) at the Very Large Telescope Interferometer (VLTI) 
is ideal for this goal with its "spectro-interferometric" capability 
and has already proven its potential to study the circumstellar 
environment of cool (and also hot) evolved stars. However, such 
studies on stellar sources with MIDI are limited to Galactic ob- 
jects up to now. Despite the large distance to the LMC (50 kpc 
is adopted in the present work), the high luminosity of WOH 
G64 and its huge infrared excess make it observable with MIDI 
(F\2fim = 7 Jy). In this paper, we present MIDI observations of 
WOH G64 and radiative transfer modeling to derive the physi- 
cal properties of its dust envelope as well as the luminosity of 
the star. 



2. Observations 

2.1. N-band visibilities 

WOH G64 was observed with MIDI in the HIGH.SENS mode 
in September and November 2005 as well as in October 2007 
with two baseline configurations using the 8.2 m Unit Telescopes 
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Fig. 1. H-R diagram with theoretical evolutionary tracks and the 
observationally derived locations of WOH G64. The black solid 
lines represent the evolutionary tracks for stars with 15, 25, 40, 
and 60 Mq with Z = 0.008 from Schaerer et al. (fT993l . The evo- 
lutionary tracks with an initial rotation velocity of 300 km s -1 
for Z = 0.008 from Meynet & Maeder (120051 ) are plotted with 
the dotted lines, if available. The red solid line represents the 
newly computed track for a 25 Mq with Z = 0.01 with the 
code of one of the authors (T. Driebe). The Hayashi limit pre- 
sented in Levesque et al. ( 120071 ) is shown with the dashed line. 
The location of WOH G64 based on r eff = 3200-3400 K and 
(5-6)xl0 5 Lq derived in the previous studies is shown with the 
box. The filled circle represents the location of WOH G64 based 
on the luminosity of 2.8 x 10 5 Lq derived in the present work 
(see Sect.|H. 



(UT2-UT3 and UT3-UT4), as summarized in Table Q] (Program 
ID: 076-D0253 and 080-D0222, P.I.: K. Ohnaka). A prism with 
a spectral resolution of A/AA ^ 30 at 10 fj.m was used to obtain 
spectrally dispersed fringes between 8 and 13 pun, and each data 
set contains typically 200 scans. A detailed description of the 
observing procedure is given in Przygodda et al. (2003), Leinert 
et al. (2004), and Chesneau et al. d2D05b . Since WOH G64 itself 
is too faint in the V band (V <: 17) for the Multi-Application 
Curvature sensing Adaptive Optics (MACAO) system, we per- 
formed "off-target Coude guiding" using the nearby star 0216- 
0046676 (USNO-B 1.0 catalog) with V= 16, which is 1076away 
from WOH G64. In addition to MACAO, the InfraRed Image 
Sensor (IRIS) operating in the K band was used to keep the 
alignment of two beams as correct as possible during interfer- 
ometric observations. 

For data reduction, we use the MIA+EWS package ver.l.fQ 
(Leinert et al. [20 04 ; Jaffe 120041 ) . Since the V magnitude of the 
guide star is close to the limit of the MACAO system on UTs 
(V = 17), our MIDI observations can be severely affected by at- 
mospheric conditions such as seeing, wind speed, and coherence 
time. In order to assess the data quality, we check the acquisition 

1 Available at http://www.strw.leidenuniv.nl/" nevec/MIDI 
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Table 1. Summary of the MIDI observations of WOH G64: night, time of observation (Coordinated Universal Time=UTC), tele- 
scope configuration (Tel.), projected baseline length B p , position angle of the projected baseline on the sky (P.A.), and seeing. The 
data sets #1 and #3 were used only for extracting the iV-band spectra. 



# 


Night 


tabs 


Tel. 


B p 


P.A. 


Seeing 


Remarks 






(UTC) 




(m) 


O 






1 


2005 Sep. 18 


08:18:03 


UT3-UT4 


57.9 


85.7 


175 


Spectrum only 


2 


2005 Nov. 12 


08:22:17 


UT3-UT4 


62.4 


136.0 


07 5-07 7 




3 


2005 Nov. 13 


08:23:15 


UT2-UT3 


33.9 


77.7 


ff.' 6 


Spectrum only 


4 


2007 Oct. 22 


05:30:49 


UT3-UT4 


56.5 


77.3 


170-175 




5 


2007 Oct. 22 


06:36:17 


UT3-UT4 


59.0 


92.7 


170-175 




6 


2007 Oct. 22 


07:41:05 


UT3-UT4 


60.8 


107.5 


170-175 





0.5 



,=62. 4m, P.A.= 136.0 
P.A. = 77.3° 



,p 



» f4fl =56. 5m, 
♦ #5 _B P = 59.0m, P.A. = 92.7° 
□ #6 B P =60.8m, P.A.= 107.5° 
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Fig. 2. a: iV-band visibilities of WOH G64 derived from the data 
sets #2, #4, #5, and #6. The errors are dominated by the system- 
atic errors as described in Sect. [2] b: Uniform-disk (UD) diame- 
ter obtained by fitting the observed visibilities. 



images of each data set, which are taken before the interfero- 
metric observations, as well as the histogram of the fringe power 
of scans at each spectral channel. This fringe histogram appears 
very broad and/or asymmetric when the image quality is poor or 
the overlap of two beams is not perfect, and therefore, it can be 
used for diagnosing the data quality. We find out that the data 
sets #1 and #3 are affected by poor image quality and/or imper- 
fect beam overlap, and therefore, we use these data sets only for 
extracting the A^-band spectra. On the other hand, the other four 
data sets do not show a hint of such problems, and we derive 
visibilities and spectra from these data sets. 

We derive the interferometer transfer function at each spec- 
tral channel between 8 and 13 fim by observing calibrators listed 
in Table [2] We use the mean of the transfer function values de- 



rived from the calibrators observed on the same night as WOH 
G64, and the errors of the calibrated visibilities are first esti- 
mated from the standard deviation of the transfer function values 
at each wavelength. While the errors estimated in this manner 
are quite small (lcr as 5%), there may still be systematic errors 
discussed above, even if the data sets seem not to be severely 
affected by the overlap problem and/or poor image quality. In 
order to account for such systematic errors, we assume a mini- 
mum relative error of 15% as adopted by Preibisch et al. (2006). 

Figure [2^ shows the A^-band visibilities of WOH G64 de- 
rived from the data sets #2, #4, #5, and #6 using the EWS pack- 
age. Since the results obtained with the MIA and EWS packages 
show good agreement, we only present the result derived with 
the EWS package in the discussion below. The object is spa- 
tially resolved by our MIDI observations. The observed A^-band 
visibilities are approximately constant or slightly decrease from 
8 to ~10 jt/m, while they increase gradually longward of 10 yum. 

Qualitatively, the observed wavelength dependence of the 
visibilities can be explained by the silicate absorption feature at 
10/im (see Fig. 0). In an optically thick dust envelope, photons in 
the silicate feature originate farther out from the star compared 
to the continuum, leading to a decrease in visibility inside the 
silicate feature. The four visibilities were obtained at approxi- 
mately the same baseline length (57-62 m) but at position angles 
differing by 59°, but there is no remarkable difference among the 
visibilities within the errors. Therefore, the object shows no sig- 
nificant deviation from centrosymmetry over this range of posi- 
tion angle. The uniform-disk diameter derived from the visibil- 
ities, which is plotted in Fig. [2j>, shows a steep increase from 
8 to 10 fim (from 18 to 26 mas) with a flat portion longward of 
10 fim. However, we stress that the uniform-disk fitting is merely 
for obtaining some characteristic angular size of the object. 

2.2. N-band spectra 

We extract the absolutely calibrated A^-band spectra of WOH 
G64 from our MIDI data with the procedure described in Ohnaka 
et al. (120071 1. The calibrators observed at airmasses similar to 
WOH G64 are used for the spectrophotometric calibration and 
marked with f in Table [2] The averaged spectra are calculated 
separately from the data taken 2005 and 2007 with three (WOH 
G64, calibrator) pairs each, and the errors are estimated from the 
standard deviation of the spectra derived from these pairs. The 
extracted A^-band spectra, which are plotted in Fig. [3] are dom- 
inated by the self-absorption of amorphous silicate at 10 jum. 
The MIDI spectrum obtained in 2007 shows flux levels slightly 
lower than that observed in 2005. However, given the errors 
of the spectrophotometric calibration of MIDI data, the differ- 
ence is marginal, and we cannot draw a definitive conclusion 
about the temporal variation of the A^-band flux between 2005 
and 2007. Also plotted are the mid-infrared spectra obtained 
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with the Spitzer Space Telescope (Werner et al. 2004) and the 
Infrared Space Observatory (ISO). The Spitzer spectrum was 
obtained on 2005 January 12 with the InfraRed Spectromete]0 
(IRS, Houck et al. 2004) in the Short Low (SL) and Long Low 
(LL) modes covering from 5.2 to 38 /im with spectral resolution 
of 64-128 (Program ID: P03426, P.I.: J. H. Kastn er) and has 
recently been published by Buchanan et al. (12006b . We down- 
loaded the Basic Calibration Data (BCD) processed with the 
S15-3 pipeline from the Spitzer data archive and extracted the 
spectrum using SMART v.6.2.4 (Higdon et al. 120041 1. The sky 
background was estimated by means of "local sky" using the 
pixels outside the extraction aperture, and the spectra extracted 
at two nodding positions are co-added. The extracted spectrum 
is in very good agreement with that published in Buchanan et al. 
d2006T >. The ISO spectrum obtained on 1996 May 26 with PHOT- 
S (2.5-5 jum and 6-12 fim, A/AA * 90), which was published by 
Trams et al. (1999), was downloaded from the ISO data archive. 
We also plot (spectro)photometric data available in the literature 
(Elias et al. [19861 Roche et al. [19931 MSX, Egan et al. [20031 
IRAS Point Source Catalog). 

The MIDI spectrum obtained in 2005 is in very good agree- 
ment with the Spitzer/IRS spectrum, which was taken relatively 
close in time to our MIDI data compared to the other data. 
We also note that the MIDI spectrum was obtained with a slit 
of 075x2", while the Spitzer/IRS uses much larger slits (376 
x57" and 1075 xl68" in the SL and LL modes, respectively). 
Therefore, the agreement between the MIDI and Spitzer/IRS 
spectra means the absence of extended emission in the N band 
(whether physically associated to WOH G64 or background 
emission). On the other hand, there are some differences in the 
flux level as well as the depth of the 10 yum feature between the 
MIDI spectra and the other data. This can be due to either differ- 
ences in the apertures used for these observations or an intrinsic 
temporal variation. Given the above good agreement between 
the MIDI and Spitzer data, the differences between the MIDI 
and the other data are more likely to be (long-term) temporal 
variations of the mid-infrared flux. In fact, Roche et al. (1993) 
report that WOH G64 shows noticeable temporal variations of 
the flux level as well as the depth of the 10 /urn absorption over 
~5 years, despite the small variability amplitude in the K band 
(AK * 0.3, Wood et al. [T992l Whitelock et al. l2003l . However, 
since the MIDI spectra taken in 2005 and 2007 do not show a 
significant temporal variation, we merge the data (Af-band vis- 
ibilities as well as spectra) taken at the two epochs for our ra- 
diative transfer modeling. In addition to our MIDI spectra, we 
also include the Spitzer/IRS spectrum in our modeling, because 
it was obtained relatively close in time to our MIDI observations. 

2.3. Spectral energy distribution 

Together with the iV-band visibilities and mid-infrared MIDI and 
Spitzer/IRS spectra, we use photometric data from the optical 
to the far-infrared as observational constraints in our radiative 
transfer modeling. We collect photometric data on WOH G64 
available in the l iteratu re (Elias et al. 119861 Wood e t al. 119921 
Whitelock et al. [20031 2MASS, Cutri et al. [200 31 MACHO 
Alcock et al . I200T1 ASAS, Pojmanski et al al. [20021 [20031 
120041 l2"005al b: IRAS Point Source Catalog; SAGE, Meixn er et 
al. I20061 >. The long-term monitoring of Wood et al. ( 1992) and 
Whitelock et al. d20031 l shows that the variability of WOH G64 



2 The IRS was a collaborative venture between Cornell University 
and Ball Aerospace Corporation funded by NASA through the Jet 
Propulsion Laboratory and Ames Research Center. 



Table 2. List of calibrators used in the present work, to- 
gether with their spectral type, 12 fim flux (F12), uniform- 
disk diameter (c/ud) and the date as well as the time stamp 
(fobs) of the MIDI observations. The uniform-disk diam- 
eters were taken from the CalVin list available at ESO 
(http://www.eso.org/observing/etc/). The stars used for spec- 
tropho tome trie calibration are marked with t- 



Calibrator 


Sp. 




d\jD 


Date 


fobs 




Type 


(Jy) 


(mas) 




(UTC) 


2005 


HD33042 


K5III 


12.3 


2.79 ±0.15 


Sep. 18 


07:37:50t 










Nov. 12 


08:46:31 










Nov. 13 


08:47:36t 


HD37160 


KOIIIb 


9.4 


2.08 ± 0.20 


Sep. 18 


08:56:35t 


HD36673 


FOIb 


8.1 


1.51+0.27 


Nov. 12 


05:33:39 










Nov. 12 


07:10:25 


HD36079 


G5II 


20.3 


2.97 ±0.16 


Nov. 12 


06:24:46 


2007 


HD26967 


K2III 


13.5 


2.59 ±0.15 


Oct. 22 


04:56:51 


HD55865 


K0III 


12.4 


2.49 ±0.13 


Oct. 22 


06:08:55 










Oct. 22 


07:05:45t 



T 




O Roche et al- (1993) 

□ Elias et al. (1986) 

A MSX 

X IRAS 



7 8 9 10 11 12 13 14 

X (/xm) 

Fig. 3. N-band spectra of WOH G64 observed with MIDI, plot- 
ted together with the Spitzer/IRS and ISO/PHOT-S spectra as 
well as (spectro)photometric data available in the literature. 
The region severely affected by the ozone band (9.6-10 fim) is 
marked by the vertical lines. 



is as small as ~0.3 in the near-infrared. The infrared spectra ob- 
tained with the PHOT-S and the Short Wavelength Spectrometer 
(SWS) onboard the ISO (Vandenbussche et al. 2002} are also 
included. Since the IRAS data only give upper limits on the 60 
and 100 fim fluxes, we derive the 70 fim flux from the imag- 
ing data obtained with the Multiband Imaging Photometer for 
Spitzer (MIPS, Rieke et al. |2004l l on 2005 November 8 (Program 
ID: P20203, PL: M. Meixner). The 70 /im flux is extracted from 
the post-BCD data with an aperture radius of 35", with the back- 
ground emission estimated in an annulus between 40" and 60". 
The resulting 70 //m flux is 1.7 Jy. The interstellar extinction to- 
ward WOH G64 is uncertain. Van Loon et al. ( 1997) argue that 
while the foreground reddening toward the LMC is Ay « 0. 1 8 on 
average, there may be considerable (but patchy) extinction inside 
the LMC. However, since no systematic study on the reddening 
inside the LMC is available, we neglect the interstellar extinc- 
tion toward WOH G64 in the present work. The whole SED is 
plotted in Fig. [4^, which reveals the huge infrared excess toward 
WOH G64. 
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Fig. 4. a: SED of WOH G64 from the optical to the far-infrared. 
The photometric data are taken from the literature and vari- 
ous catalogs as given in the panel. The data from Wood et al. 
( 119921 ) and Whitelock et al. (2003 ) are the time-averaged val- 
ues. Only the 5.8 /mi flux measured with the InfraRed Array 
Camera (IRAQ of Spitzer is currently available in the SAGE 
database. The IRAS 60 //m and 100 pm fluxes represent up- 
per limits, b: The H2O absorption features in the region be- 
tween 1 and 8 pm. Solid line (between 2.7 and 4 pm): ISO/SWS. 
Crosses: ISO/PHOT-S. Filled circle s: Spit zer/IRS. Open circles: 
Photometric data from Wood et al. (|1992| l. 

Figure [4J5 shows an enlarged view of the region between 1 
and 8 pm. The ISO/SWS and ISO/PHOT-S spectra reveal re- 
markable absorption features between 2.4 and 4 /mi, which can 
be attributed to the H2O v\ and V3 fundamental bands. As pre- 
sented in Sect. [3] the flux contribution of the attenuated star is the 
more dominant compared to the dust thermal emission at wave- 
lengths below ~3 pm. Therefore, the H2O absorption feature 
centered at 2.7 pm is likely to originate in the photosphere and/or 
in the so-called "warm molecular layers" (or MOLsphere) close 
to the star, whose presence is confirmed in Galactic RSGs (e.g., 
Tsuji 2000] [20061 Ohnaka |2004t Verhoelst et al. [20061 Perrin 
et al. 120071 ). On the other hand, the Spitzer/IRS spectrum (and 
also the ISO/PHOT-S spectrum, although a little noisy) reveals 
the 6 pm H2O feature due to the V2 fundamentals, which is only 
marginally discernible in Fig. la of Buchanan et al. (2006) be- 
cause of the large ordinate scale of their plot. Since the flux is 
dominated by the dust thermal emission at these wavelengths, 
the 6 pm absorption feature is mostly caused by the H2O gas 
present in the circumstellar envelope, and the H2O absorption is 
seen with the hot dust thermal emission as the background. Such 



circumstellar H2O V2 absorption features, albeit much weaker, 
are also identified in the Galactic red supergiant NML Cyg 
(Justtanont et al. |19961 ). 

3. Modeling of the dusty environment of WOH G64 

Since van Loon et al. ( 119991 ) have already shown that the 
SED of WOH G64 cannot be explained by a spherical dust 
shell model, we attempt to explain the observed SED and N- 
band visibilities with the second simplest geometry: axisymmet- 
ric dust torus (or disk) model. We use our Monte Carlo code 
(mcsimjnpi) for two-dimensional modeling of the dusty envi- 
ronment of WOH G64. The outputs of mcsimjnpi are dust tem- 
perature and monochromatic mean intensity in each cell of the 
dust envelope. Using these outputs, a model image at a given 
wavelength viewed from an arbitrary angle is computed with our 
ray-trace code. The details of our codes are described in Ohnaka 
et al. (|2006l ). 

The central star is assumed to be a blackbody of effective 
temperature T e ff . The spectral type of WOH G64 is determined 
as M7.5 from optical spectroscopy by Elias ( 1986 ) and van Loon 
et al. ( 120051 ), although the latter authors note some ambiguity in 
the spectral classification of WOH G64 and suggest that an ear- 
lier spectral type of M5 is also possible. We convert these spec- 
tral types to effective temperatures using the temperature scale 
derived by Levesque et al. ((2006) based on RSGs in the LMC. 
Since their effective temperature scale does not extend beyond 
M4.5, we linearly extrapolate the temperature scale between M0 
and M4.5 down to M7.5, which results in r e (j = 3230 K. The ear- 
lier spectral type of M5 corresponds to T e ff = 3400 K. We adopt 
r e ff = 3200 K in our modeling, but the adoption of 3400 K turns 
out not to noticeably affect the result of the modeling presented 
below. 

The torus model considered here is similar to that adopted by 
Ohnaka et al. (2006), which is basically a sphere with a biconical 
cavity. We introduce two differences compared to this previous 
model: 1) the density in the cavity is lower than in the torus but 
not exactly zero, and 2) the density drop from the torus to the 
cavity (along the latitudinal direction) is smoothed as described 
by 

p(r , 0) oc X ( e(|9| _ | /£+1 + e _ M _e )/e+1 ) , 

where r and 8 are the radial distance from the central star and 
the latitudinal angle measured from the equatorial plane, respec- 
tively, while r m and are the inner boundary radius and the 
half-opening angle of the torus, respectively (see the inset of 
Fig. [5^). e is a constant representing the smoothing of the den- 
sity drop from the torus to the bipolar cavity. In the models dis- 
cussed below, we adopt e = 0.03 (for 8 and measured in ra- 
dian). / represents the ratio of the density in the cavity to that in 
the torus. The radial density distribution is assumed to be pro- 
portional to r~ 2 (i.e., p = 2), which corresponds to a steady 
mass outflow with a constant velocity. The free parameters in 
our models are r m , 0, -ry (the radial optical depth at 0.55 pm 
in the equatorial plane), and /. In addition to these parameters 
describing the torus, it is necessary to specify its orientation in 
the plane of the sky. This is equivalent to specifying the angle 
between the torus symmetry axis projected onto the plane of the 
sky and the projected baseline vector of our MIDI observations 
(see also Fig. [6), and this angle is used as an additional free pa- 
rameter. The torus is assumed to consist of amorphous silicate 
grains, because the MIDI, Spitzer/IRS, and ISO spectra show no 
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obvious dust features other than the self-absorption due to amor- 
phous silicate. The grain size distribution is proportional to aT q 
between the minimum size (a m i„) and maximum size (a max ), as 
considered by Mathis et al. dl977l l. We tentatively assume a m i n 
= 0.005 fim, a max = 0.1 /mi, and q = 3.5, and examine the ef- 
fects of these parameters later. The mean of the absorption and 
scattering coefficients weighted with the above size distribution 
is computed from the complex refractive index of the "warm sil- 
icate" of Ossenkopf et al. ( 1 19921 ) using the Mie code of Bohren 
& Huffman ( 119831 1. For simplicity, isotropic scattering is consid- 
ered for our modeling. 

The range of the free parameters varied in our modeling is 
as follows: r m (R+) = 5 ... 25 with Ar in = 5, = 20° ... 90° 
with A0 = 10° (0 = 90° corresponds to spherical shell mod- 
els), t v = 10 ... 50 with At v = 5, and / = 0.0 ... 0.5 with A/ = 
0. 1 . The outer boundary radius (r out ) is fixed to 250xr; n (see dis- 
cussion below). SEDs are computed for different inclination an- 
gles (measured from pole-on as illustrated in Fig.|5^) of i = 0°... 
90° with Ai =10° for all models, and we first compare these 
model SEDs with the observed one. This comparison between 
the model and observed SEDs allows us to derive the bolometric 
luminosity of the central star, because the distance to the LMC is 
known (50 kpc). The two-dimensional visibilities are then com- 
puted from the Fourier transform of model iV-band images. Note 
that the stellar radius derived from the central star's luminosity 
and effective temperature is used to scale the model visibilities. 

Figure [5] shows a comparison of the best-fit model with the 
observed SED and A^-band visibilities, and the parameters of the 
model are given in Table [3] As shown in the inset of Fig. [5^, the 
best-fit model is characterized by a geometrically and optically 
thick torus viewed close to pole-on. The whole SED including 
the 10 fim silicate feature is reasonably reproduced. The torus 
is optically thick in the equatorial direction (jy = 30 and Tio^m 
= 7), which gives rise to the huge mid- and far-infrared excess 
as observed. However, since the object is seen through the cav- 
ity with lower density (tv = 9 and Tw^m - 2.1), the central star 
is not entirely obscured, and the flux in the near-infrared and in 
the optical as well as the observed silicate absorption is reason- 
ably reproduced. The model predicts the near-infrared (H and 
K bands) flux to be somewhat lower than the observation. This 
is because the spectrum of the central star is approximated with 
that of a blackbody. The spectrum of photospheric models shows 
a flux excess at these wavelengths compared to the blackbody of 
the same effective temperature, because the opacity of negative 
hydrogen, which is the major continuous opacity source in the 
atmosphere of cool stars, has a minimum at ~1.6 fim (e.g., Tsuji 
119781 Plez et al. 11992) . The discrepancy in the optical flux be- 
tween the model and the observation can be explained by molec- 
ular absorption features (mainly TiO) which are not included in 
the blackbody approximation. For comparison, we show a spher- 
ical shell model with the same parameters as the best-fit torus 
model, but with -ry = 20 (and a higher luminosity of 3.8 x 10 5 Lq 
than the 2.8 x 10 5 Lq of the torus model discussed below). This 
spherical shell model can reproduce the mid- and far-infrared 
excess as well as the 10 fan silicate feature, but the predicted 
near-infrared and visible fluxes are too low to explain the ob- 
served data, because the optical depth along the line of sight is 
much higher. 

The model visibilities are also in agreement with the MIDI 
observations within the errors of the measurements. This pole- 
on model appears nearly centrosymmetric as seen in the A^-band 
images shown in Fig. [6] and this can explain the absence of a no- 
ticeable position-angle dependence in the observed visibilities. 
The orientation of the torus in the plane of the sky cannot be 
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Fig. 5. Best-fit dust torus model for WOH G64. As shown in 
the inset of the panel a, the torus is geometrically and optically 
thick and viewed relatively close to pole-on (see Table[3]for the 
parameters), a: Comparison of the observed and model SEDs. 
The SED of the torus model is plotted by the thick solid line. The 
symbols as well as the blue, green, and red lines represent the ob- 
served SED (see Fig.|4^ for the references of the symbols). The 
dotted line represents the stellar flux contribution of the torus 
model. A spherical shell model is plotted with the dashed-dotted 
line. This spherical model has ry = 20, and the other parame- 
ters are the same as the best-fit torus model, b: Comparison of 
the mid-infrared spectra, c: The observed A^-band visibilities are 
plotted by the filled circles (data set #2), filled triangles (#4), 
filled diamonds (#5), and open squares (#6). The solid, dashed, 
and dotted lines (almost entirely overlapping with one another) 
represent the model visibilities calculated for different orienta- 
tions of the torus symmetry axis (solid line: torus axis perpen- 
dicular to the baseline vector, dashed line: 45°, dotted line: par- 
allel), as illustrated in Fig. [6] These model visibilities show no 
notable difference, because the model is viewed close to pole-on. 
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Fig. 6. iV-band images predicted by the best-fit torus model for WOH G64 shown in Fig. The color scale is normalized with 
the maximum intensity (excluding the star) and the minimum intensity set to 1 % of the maximum intensity. The inner rim of the 
torus can be seen as the bright (red) circle on the 8 /mi image. The orientations of the baseline vector used for the calculation of the 
visibilities in Fig.[5j; are shown by the solid, dashed, and dotted lines. 



well constrained because the object is seen close to pole-on, but 
obviously, it is not a crucial parameter in such a case and does 
not affect the determination of the other parameters. We note 
that the inner boundary radius is weakly constrained, 10-20 R+, 
because the four visibilities currently available were taken at ap- 
proximately the same baseline length and cannot put tight con- 
straints on the radial density distribution. This hampers us from 
discussing possible metallicity effects on dust formation in LMC 
RSGs at the moment. 

We study the effects of the outer boundary radius by vary- 
ing r out in the best-fit model. While the adoption of larger outer 
boundary radii (r out = 500 r m and 1000 r m ) does not result in 
a notable change in the SED below 70 yum and A^-band visi- 
bilities, the models with the outer boundary radii smaller than 
250 r m predict the 70 fim flux to be too low compared to the 
Spitzer/MIPS data. Therefore, we can only set a lower limit of 
250xr; n on the outer boundary radius. We also check the effects 
of the grain size distribution by changing the minimum and max- 
imum grain sizes as well as the power-law index of the size dis- 
tribution for the best-fit model. We find the upper limit for the 
minimum grain size to be 0.05 yum, while the lower limit can- 
not be well constrained. On the other hand, the upper and lower 
limits for the maximum grain size are found to be ~0.15 /im 
and ~0.01 /urn, respectively. Models with power-law index val- 
ues between 1.5 and 4.5 (3.5 is adopted in the above modeling) 
can reproduce the observational data, and thus, this parameter is 
not tightly constrained by the present modeling. 

4. Discussion 

Our two-dimensional modeling of the MIDI data and SED indi- 
cates that WOH G64 has an axisymmetric, optically thick dust 
torus. The presence of such axisymmetric structures is suggested 
toward Galactic RSGs as well. In particular, the well-studied 
dusty RSG, NML Cyg, which can be regarded as a Galactic 
counterpart of WOH G64 (van Loon et al. 119981 1, is also sus- 
pected to have bipolar outflows based on radio interferometry 
observations (Richards et al. 1996). Therefore, the occurrence 
of bipolar outflows and/or thick tori seems to be a common phe- 



Table 3. The parameters of the circumstellar dust torus and 
the central star of WOH G64 derived from our radiative trans- 
fer modeling. 



Parameter 



Value 



Torus half-opening angle (0) 
Torus inner boundary radius (ri n ) 
(with 1 R+ = 8 AU) 
Inner boundary dust temperature 
Torus outer boundary radius (r out ) 
Radial optical depth at 0.55 fim (t v ) 
Optical depth along the line of sight 
Radial density distribution 
Inclination angle (0 
Grain size distribution 
Maximum grain size (a max ) 
Minimum grain size (a min ) 
Total envelope mass 

Luminosity 
Stellar radius 
Effective temperature 



60° ±10° 
15 ±5tf* 
(120 ±40 AU) 
880!|* K 
> 250 X r m 
30 ±5 

9±2(/ = 0.3+0.1) 
r~ 2 (fixed) 

20° (0°- 30° with i < 90° - 0) 

-3 S 

oc a 

0.1 urn (0.01-0.15 /mi) 
0.005 yum (<0.05 /im) 
3-9 M Q 

(gas-to-dust ratio = 200-500) 
(2.8 ± 0.3) x 10 5 Lq 
1730 R e 
3200 K (fixed) 



nomenon among RSGs in the Galaxy as well as in the LMC, 
although the origin of the axisymmetric structures is still un- 
clear. We note that the pole-on dusty torus model for WOH G64 
is qualitatively consistent with the asymmetric OH maser profile 
with the blue peak much stronger than the red peak. Marshall 
et al. (|2004| l propose that the amplification of the stellar light 
is responsible for the blue-red-asymmetry, and this process can 
be particularly effective, when the star is seen through the less 
dense cavity as in the case of WOH G64. 

The luminosity of the central star derived from our radia- 
tive transfer modeling, 2.8 x 10 5 Lq (Mb i = -8.8), is by a fac- 
tor of ~2 lower than the (5-6) xlO 5 Lq previously estimated by 
Elias et al. (1986]! and van Loon et al. (2005) assuming a spheri- 
cal shell. Our modeling suggests that the object has an optically 
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thick torus and is seen through the cavity with lower density. In 
such a case, more radiation escapes toward the cavity (i.e., to- 
ward the observer) than toward the optically thick torus, and the 
bolometric luminosity of the central star is overestimated, when 
it is derived from the observed flux assuming isotropic radiation. 
The location of WOH G64 based on the new, lower luminosity 
is plotted in Fig. Q] The newly derived luminosity corresponds 
to that of an RSG with an initial mass of ~25 Mq, and the loca- 
tion of WOH G64 is in much better agreement with the theoret- 
ical evolutionary tracks, although WOH G64 is still somewhat 
too cool compared to the theoretical models. In fact, compari- 
son with the Hayashi limit presented in Levesque et al. (2007 ) 
reveals that WOH G64 is to the right of the Hayashi limit if T e fi 
is 3200 K. Even if we adopt T e ff = 3400 K corresponding to the 
earlier spectral type of M5 suggested by van Loon et al. (2005), 
it is close to, but still slightly beyond the Hayashi limit. The ex- 
treme effective temperature of WOH G64 suggests that it may be 
experiencing a very unstable phase accompanied by heavy mass 
loss, as proposed for some RSGs in the LMC and SMC located 
in the forbidden zone beyond the Hayashi limit (Levesque et al. 
2007 ). However, it should also be kept in mind that the determi- 
nation of effective temperature from the optical spectrum may 
be affected by the extinction due to the circumstellar dust enve- 
lope and/or the contribution of scattered light by dust grains, as 
suggested from the redshifted photospheric absorption lines with 
respect to the SiO and H2O masers (van Loon et al. 1998, 2001 ). 
Also, the Hayashi limit is a complex function of stellar mass and 
metallicity as described in Cox & Giuli ( 119681 ). Therefore, it is 
still possible that WOH G64 is just on the Hayashi limit with an 
effective temperature somewhat higher than the above estimates. 

WOH G64 indeed appears to have already shed a significant 
fraction of its initial mass. The dust mass in the best-fit torus 
model is ~ 1 .7 x 10~ 2 Mq with a grain bulk density of 3 g cm" 3 , 
which translates into a total envelope mass of ~3-9 Mq with 
gas-to-dust ratios of 200-500. These values already represent 
a noticeable fraction of the initial total mass of ~25 Mq es- 
timated for WOH G64. Moreover, as discussed in Sect. [?] the 
outer boundary radius of 250xn n is a lower limit, and the enve- 
lope mass becomes even higher for larger outer boundary radii. 
Interestingly, although WOH G64 is cooler than the evolution- 
ary tracks for a 25 Mq star shown in Fig. [1] the evolutionary 
calculation from Schaerer et al. {T993 ) as well as the new calcu- 
lation performed in this work predicts that a 25 Mq star has lost 
2 and 7.5 Mq at the RSG stage (an age of 7.8 Myr), respectively. 
Therefore, the derived envelope mass of WOH G64 is in rough 
agreement with current stellar evolution theory. 

5. Concluding remarks 

We have spatially resolved the circumstellar environment of the 
dusty red supergiant WOH G64 in the LMC, which is the first 
MIDI observation to resolve an individual stellar source in an ex- 
tragalactic system. The observed A/-band uniform-disk diameter 
is found to increase from 18 mas at 8 fim to 26 mas at 13 fim. The 
visibilities obtained at approximately the same baseline length 
but at four position angles differing by ~60° show no remarkable 
difference, indicating that the object appears nearly centrosym- 
metric. Our two-dimensional radiative transfer modeling shows 
that these visibilities and the observed SED can be explained by 
an optically and geometrically thick torus viewed close to pole- 
on. The luminosity of the central star derived by our modeling 
is by a factor of ~2 lower than the previous estimates based on 
spherical shell models. The new, lower luminosity brings WOH 
G64 in much better agreement with the theoretical evolutionary 



tracks, which suggests an initial mass of ~25 Mq for WOH G64. 
However, its location on the H-R diagram is very near or beyond 
the Hayashi limit. This implies that WOH G64 may be in an un- 
stable stage at the moment, experiencing remarkable mass loss. 
The total envelope mass estimated from our model is of the or- 
der of 3-9 Mq, representing a significant fraction of its initial 
mass. This envelope mass is comparable to what is predicted by 
current stellar evolutionary calculations. 

We also identified the H2O absorption features at 2.7 jum 
and 6 fim in the spectra obtained with ISO and Spitzer, respec- 
tively. The 2.7 fim feature is likely to originate in the photosphere 
and/or MOLsphere. On the other hand, the 6 fim feature can be 
interpreted as the absorption due to cold circumstellar H2O gas 
seen with the hot dust thermal emission as background. 

Obviously, follow-up MIDI observations are necessary for 
putting more constraints on the properties of the torus. In par- 
ticular, MIDI observations with shorter and longer baselines are 
indispensable for obtaining tighter constraints on the torus in- 
ner boundary radius, which is crucial for better understanding 
possible metallicity effects on the dust formation and mass loss 
in RSGs. MIDI measurements over even a wider position an- 
gle range are also important for probing the possible presence 
of more complex structures. The H2O gas in the photosphere 
and in the MOLsphere can be studied in more detail using the 
near-infrared VLTI instrument AMBER (Petrov et al. 2007). For 
example, AMBER observations of the H2O bands in the K band 
would be useful for constraining the geometrical extension of 
the H2O layers close to the star. 
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